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ABSTRACT 

Wc study the dynamics of galactic disk formation and evolution in 'realistic' A cold 
dark matter haloes with idealized baryonic initial conditions. We add rotating spheres 
of hot gas at z = 1.3 to two fully cosmological dark-matter-only halo (re)simulations. 
The gas cools according to an artificial and adjustable cooling function to form a rota- 
tionally supported galaxy. The simulations evolve in the full cosmological context until 
z=0. We vary the angular momentum and density profiles of the initial gas sphere, 
the cooling time and the orientation of the angular momentum vector to study the 
effects on the formation and evolution of the disk. The final disks show exponential 
radial and (double)-exponential vertical stellar density profiles and stellar velocity dis- 
persions that increase with age of the stars, as in real disk galaxies. The slower the 
cooling/accretion processes, the higher the kinematic disk-to-bulge (D/B) ratio of the 
resulting system. We find that the initial orientation of the baryonic angular momen- 
tum with respect to the halo has a major effect on the resulting D/B. The most stable 
systems result from orientations parallel to the halo minor axis. Despite the spherical 
and coherently rotating initial gas distribution, the orientation of the central disk and 
of the outer gas components and the relative angle between the components can all 
change by more than 90 degrees over several billion years. Initial orientations perpen- 
dicular to the major axis tend to align with the minor axis during their evolution, but 
the sign of the spin can have a strong effect. Disks can form from initial conditions 
oriented parallel to the major axis, but there is often strong misalignment between 
inner and outer material. The more the orientation of the baryonic angular momentum 
changes during the evolution, the lower the final D/B. The behaviour varies strongly 
from halo to halo. Even our very simple initial conditions can lead to strong bars, 
dominant bulges, massive, misaligned rings and counter-rotating components. We dis- 
cuss how our results may relate to the failure or success of fully cosmological disk 
formation simulations. 
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1 INTRODUCTION 

The majority of galaxies in the local universe with masses 
similar to that of the Milky Way are disk-dominated (e.g. 
iDelgado-Serrano et al.l |2010|). In the standard paradigm, 
these disks formed through cooling and condensat i on of 
gas within cold dark m atter haloes IjWhite &i ReesI 119781 : 
iFall fc EfstathioulllQSd ). Disk size is a consequence of the 
angular momentum of the gas, which had previously been 
acqu ired through t idal torques from neighboring structures 
(e.g. iPeeblell 19691 : IWhitelll984l ). 

Numerical hydrodynamical simulations of this for- 
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mation scenario have been carried out in great number 
(e.g. Navarro & White 1994; Navarro & Stoinmctz 1993; 
Abadi et al.ll2003l : iGovernato et al.ll2007l: IScannapieco et aU 



2009|). Despite recent progress (e .g. iGuedes et al.l I2OI1I : 



Agertz et al.ll201ll : ISales et al.ll201ll ). these simulations have 
in general suffered from a range of problems, includ- 
ing angular momentum loss leading to overly small disks 
IINavarro fc Benall99ll ) , production of overly massive galax- 
ies ( Guo et al.l I2OI0I), and too much early and too little 



late star formation IjScannapieco et al.ll2009l ). Even if ex- 
tended di sks form, th ey are often destroyed by infalling satel- 
lites (Toth fc Ostriker 1992) or the accretion o f misaligned 
gas IjScannapieco et al.ll2009l : ISales et al.ll201lh . Moreover, 
different numerical schemes can yield very different galax- 
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ies for the same i nitial conditions (iScannapiec o et al.ll201ll : 
iKeres et al.1 I2OIII ). In particular, the angular momentum 
content and thus the structure of the disks depends s trongly 
on the numerical method applied (JTorrev et al.ll201H ) . Over- 
all, the galaxy population predicted by cosmological gas dy- 
namical simulations disagrees in major ways with observa- 
tions. Since detailed population properties are reproduced 
quite well in the standard ACDM paradigm by simple semi- 
analytic simulations (e.g. IGuo et al.ll201H ). it has been ar- 
gued that the relevant star formation and feedback processes 
are still inadequately represented in hydrodynamical simu- 
lations. Such semi-analytic models are however not capable 
of properly capturing the complex dynamical interactions 
between gas, stars and dark matter found to be important 
in the above-mentioned simulations. 

In contrast, there is general agreement on the forma- 
tion and structure of dark matter haloes in ACDM (e.g. 
ISpringel et al.ll20(M ). and on the observational side, a de- 
tailed picture of the structure of disk galaxies has been as- 
semb led over the last decades (see lvan der Kruit fc FreemanI 
I2OIII for a recent review). In this study, we therefore study 
whether the detailed output of simulations of dark matter 
halo formation and our detailed knowledge of disk galax- 
ies can be brought into agreement if the fully cosmological 
treatment of baryon physics is replaced by idealized models. 

P revious att e mpts in this direction include the mod- 
els of IWeil et al.l ([1993) > who showed that the prevention 
of radiative gas cooling in cosmological hydrodynamical 
simulations until z = 1 allows the formation of a disk 
galaxy populatio n with realist ic angular momentum content. 
iKaufmann et al.l l|2006l . 120071 ) simulated disk formation by 
allowing a rotating gas distribution to cool inside idealized 
(mostly spherical) equilibrium haloes. This enabled them to 
examine numerical effects as well as some physical processes 
related to disk simulations. 

ACDM haloes, however, show substructure, are triax- 
ial and are continuously accreting (|Frenk et al.lll98a) . The 
impact of halo shape is still an open issue. Disk galax- 
ies in cosmological hydrodynamical si mulations are typ i- 
cally aligned with the halo minor axis IjBailin et al.ll2005h . 
Moreover, it has been argued that the interaction of the 
forming disk galaxy with the dark matter tends to make 
the haloes axisymmetri c (e.g. Berentzen fc Shlosmanl 20061 : 
iKazantzidis et aLllJOlol ). iBerentzen fc Shlosmanl l|2006l ') in- 
troduced a disk galaxy within a cosmological triaxial halo 
by gradually adding stellar disk particles according to an 
axisymmetric analytical disk model with growing mass and 
size, to the simulation. They were thus able to analyze the 
interplay between a galactic disk and a triaxial halo. A sim- 
ilar study is currently being undertaken for the same dark 
matter haloes we use here by DeBuhr et al. (in prep.). 

For our study, we perform a set of controlled numerical 
experiments using simplified prescriptions for the physics of 
gas accretion and star formation in order to gain insight into 
the dynamical processes that affect the formation of galaxy 
disks within ACDM haloes . We use fully cosmological, triax- 
ial haloes as initial conditions, but simulate disk formation 
and evolution by the cooling of a rotating gas sphere start- 
ing at redshift 2 = 1.3. Our idealized treatment allows us 
to study how the formation of a galactic disk is affected by 
the rapidity of gas cooling and by the angular momentum 
and mass of the gas and by the orientation of its angular 



momentum with respect to the principal axes of the dark 
halo. We compare our results to the direct output of fully 
cosmological hydrodynamical simulations carried out for the 
same haloes, so that we can better understand why they fail 
to produce substantial disks. 

In Section 2 we describe the setup and numerical meth- 
ods of this work. In Section 3 we describe the formation, 
evolution and structure of one particular disk model. In Sec- 
tion 4 we analyze the dependence of the disk growth on our 
model parameters. In Section 5 we discuss the influence of 
initial spin orientation and the stability of the orientation 
of the resulting disks. In Section 6 we compare our mod- 
els to the fully cosmologi cal galaxy formation simulations of 
IScannapieco et aU l|2009l ) (CS09 hereafter). Finally, in Sec- 
tion 7 we summarize and conclude. 



2 SIMULATION SETUP 

To mimic the formation of a galactic disks within fuUy- 
cosmological dark-matter-only simulations of haloes ex- 
pected to host Milky- Way type gala xies, we make use of the 
Aquarius simulations IjSpringel et al. 2008), a suite of high 
resolution zoom-in resimulations of six dark matter haloes. 
These haloes were chosen from a simulation of a cosmolog- 
ical box with a side-length of 137 Mpc, and were simulated 
from z = 127 assuming a ACDM universe with the following 
parameters: JIa = 0.75, f2m = 0.25, Q,b = 0.0 4, gg = 0.9 and 
Hq — 73 kms~^Mpc~^. For details we refer to lSpringel et al.l 
(i2008il . 

We select two haloes, named A and C, both at resolution 
level 5, which corresponds to dark matter particle masses of 
mdm « 3 X 10^ Mq and a gravitational softening length of 
fdm = 685 pc. These haloes have (dark matter only) z = 
virial masses of M200 = 1.853 and 1.793 x 10^^ M© and have 
been studied in the fully cosmological galaxy formation sim- 
ulations of CS09. We selected halo A as it has a very qui- 
escent merger history after redshift z ~ 1 and has been 
identified as a prime candidate for hosting a disk galaxy by 
semi-analytical modeling. However, CS09 did not find a sig- 
nificant disk component in their z = Q simulations. Halo C, 
in contrast, showed the highest (yet still unrealistically low) 
disk-to-bulge ratio, D jB ~ 1/4, of all 8 haloes studied in 
this work. H alo C was also studied in the Aquila Compari- 
son Project IjScannapieco et al.ll201ll ). which compared the 
results of various cosmological gas-dynamical codes. 

For our simulations w e apply the TreeSPH-code 
GADGET-3, last described in lSpringell l|2005l ). an extended 
version of which was used by CS09. 

As initial conditions for our nu mer ical experiments 
we use the simulation outputs of iSpringel et a l. (2008) at 
z — 1.3. We choose this time, as it falls after the epoch 
of major mergers. Disk galaxies have very likely been in 
place and continuously forming since before z = 1.3. For 
example, age es timates for th e solar neighbourhood exceed 
10 Gyr fsee iAumer fc Binnevi 2009 and references therein), 
making our choice appear problematic. It has been argued 
that the ages of th in disk stars are con sistent with formation 
after z ~ 1.5 (e.g. iReddv et al.ll2006l ). but it is not clear, if 
a two-component-division of the Milky Way disk into thin 
and thick components is sensible (e.g. Schonrich fc Binnevi 
[20091). Our goal is not however to model the entire forma- 
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tion and evolution of disk galaxies, but to study the ability 
of realistic dark matter haloes with evolving substructure to 
host thin disk galaxies over cosmological times. 

To simulate the formation of a disk galaxy within 
the dark matter halo, we insert a rotating sphere of gas 
into the inner halo and evolve the combined system until 
z = 0. In our models, the gas component is initially hot 
(r ~ 10® K), with an internal energy structure determined 
by the assumption of approximate hydrostat ic equilibrium 
fcf iKaufmann et al.|[200g : lAumer et al.ll2010h . The pressure 
profile is thus 



fT?\ I GMtot(r) , 
P{R) = / Paas 5 dr, 



(1) 



where Mtot{r) is the spherically averaged mass of dark mat- 
ter and gas within spherical ra dius r. The gas sphere has a 
density profile p(r) oc r^^ fe.g. iNavarro fc Whitdll993h and 
its total mass is Mtot = 5.0 — 7.5 x IO'^'^Mq, in the range 
of es timated stellar masses for the Milky Way (IMcMillanl 
uOllh and of galaxy masses expected for this halo mass 
range by abundance matching techniques (e.g. IGuo et al.l 
I2OI0I ). We choose this profile as tests of various other profiles 
yie lded undesired arti facts. An NFW-like profile, as used 
bv lAumer et al.l (|2010r ). leads to more dominant dispersion- 
dominated populations forming from the initially greater 
central gas mass, which adjusts to the triaxial potential of 
the central halo before transition to a disk-like configuration. 
For a constant density sphere, however, the central build-up 
of baryonic mass is initially too slow, delaying the transition 
from triaxial to disk-like potential relative to the baryonic 
mass build-up. 

The gas sphere is truncated at radius -Rgas = 60 — 
100 kpc. We assume zero pressure outside the gas sphere 
when determining the pressure profile according to equation 
[1] This leads to an initial expansion phase for the outermost 
gas layers due to the finite temperatures assigned to the 
SPH particles. However, this effect is negligible compared 
to the adjustment of the initially spherical gas distribution 
to the triaxial and substructured gravitational potential of 
the dark halo, which prevents the setup of idealized equi- 
librium initial conditions. The effects of this adjustment are 
discussed in Sections 3 and 4. Still, as we also show in these 
sections, this simple setup is justified by the fact that it leads 
to disk galaxy models with realistic properties. 

We use a baryonic particle mass of TTt bar ~ 10 Mp and 
a gra vitational softening ebar ~ 205 pc. IKaufmann et al.l 
(|2007. ) found that the angular momentum content and mor- 
phology of forming disks depend on resolution. Our resolu- 
tion is similar to their highest resolution, for which they 
found significantly reduced angular momentum losses. In 
section 4, we show that angular momentum losses due to 
lack of resolution are negligible in our simulations. The gas 
is set to rotate with a rotation velocity profile Kot(r), which 
is either constant or changing linearly from central to outer 
spherical radii r. These profiles yield realistic disk mass pro- 
files as is shown in the following sections. The angular mo- 
mentum vector of the gas is aligned either with one of the 
principal axes of the halo or with its angular momentum in 
order to study how this affects the final orientation of the 
disk (see section 5). 

To model the cooling of the hot gas into the centre of 
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Figure 1. The ratio of the cooling time to the free-fall time 
*cool/tff S'S a function of radius for several models with different 
initial density and cooling time parameters. AExtlC, AExt2C 
and AExt only differ by normalization of the cooling time profile 
*cool 0- AExtLoMa and AExt2C only differ in total gas mass Mgas. 
AExtLoMa and ARef only differ in initial radius of the gas sphere 
ijgas. For a full list of models and their properties see table [Tl 



the halo, we apply a simple parametrized cooling function 
with 



f-cool — Ccool.O 



Pcool.O 



(2) 



where the normalization fcooi.o determines the rapidity of 
coohng. We choose fcooi.o = 0.5 — 2.0 Gyr for pcooi.o = 
IO^'Mq kpc"^ « 4 X 10"''mH cm"^ and a = 0.56. These 
values were chosen to combine a simple dependence of the 
cooling time on the density with cooling time values which 
initially are a few times the free fall time tg at all radii, 
^cooi increasing with radius, yet the ratio icooi/iff decreasing 
with radius. This enables an initially slow transition of the 
triaxial dark-matter dominated centre of the halo to a disk- 
dominated quasi-axisymmetric system. Ratios of icooi/iff as 
a function of radius for several models with different initial 
density and cooling time parameters are shown in figure[l]to 
illustrate how these parameters affect the cooling of gas in 
our models. The initial cooling time profiles are affected by 
the implications of the non-equilibrium initial conditions as 
discussed above on dynamical timescales and by the cooling 
of gas and star formation on cooling timescales. In Sections 
3 and 4, we show that this cooling time setup produces gas 

accretion histories appropriate for our mo dels. 

A temperature floor Tfloor = 10^ K (jKaufmann et al.l 

[200y) is applied to ensure that the fo rming disk does no t 
fragment into a few massive clumps IjAumer et al.l I2OI0I ) , 
which would then enhance bulge formation. The high choice 
of Tfloor = 10^ K also prevents dense clump s from forming in 
the in falling gas via the cooling instability IjKaufmann et al.l 
l2006h . 

Such a setup clearly ignores evidence from cosmolog- 
ical hydrodynamical simulations that steady, narrow, cold 
gas streams penetrating the shock-heated atmosphere of 
massive dark matter haloes, are important in feeding gas 
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Table 1. Overview over our models 

Column 1 :halo A or C; Column 2:our model name; Column 5:the radius of the initial gas sphere; Column 4 :the mass of the initial gas 
sphere; Colum,n 5:the rotational velocity of the initial gas sphere, a single value stands for a constant Vrot{r), a — b kms~^ stands for a 
linear increase in V'rot('') from a in the centre to b at iJgas; Column 6:the number of baryonic particles; Column 7:the cooling time profile 
applied; the number stands for the approximate ratio icool/*dyn(^gas) (see text and figure[TJ; Column S:the orientation of the model. 
AM300 stands for the orientation of the halo angular momentum within 300 kpc, Major, Medium and Minor stand for the principal 
axes of the halo, '-I-' and '-' are arbitrarily chosen and differ by 180 degrees, EndARef stands for the z = orientation of model ARef, 
45 stands for an orientation halfing the angle between Major and Minor, Cosmo stands for the orientation of galaxy angular momentum 
in the cosmological runs of CS09, CArb is an arbitrary orientation in halo C, EndCExt-|-Min stands for the z = orientation of model 
CExt-fMin. 



onto forming galaxies fe.g. lPekel et al.ll2009! ). However, cold 
streams are onl y expected to be the domin ant accretion 
mode at z > 2 (|Faucher-Giguere et al.ll201ll ) and even at 
these redshift s, there is little dir ect observational evidence 
for them (e.g. ISteidel et al.ll2010l ). The majority of the stel- 
Isj ma;SS in the M ilky Way disk was formed after 2 ~ 2 
IjAumer fc Binnevli2009 ) and it may be argued that the un- 
derlying continuous star formation since z ~ 1 is driven by 
cooling of hot coronal gas in the wakes of galactic foun- 
tain clouds, which transf ers gas from the virial- temperature 
corona to the disk (e.g. iMarinacci et al.ll20ld ). Our focus 
here is not on simulating a realistic assembly process, but 



rather on using a simplified scheme, which allows important 
aspects of the underlying dynamics to be clarified. 

In our simulations, once the gas has reached the tem- 
perature floor and crossed a density threshold po.sfr, it can 
form stars with 



f-sfr — ^s 



PO.sfr 



(3) 



We choose po.sfr to correspond to a hydrogen number den- 
sity nn, thresh = 0.1 cui""^ and we use /? = 0.5 and fsfr.o ~ 
10 X idyn(po,sfr), as ofteu applied in simulations of this res- 
olution (CS09). Gas particles eligible for star formation are 
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stochastically turned into colli sionless particles of the same 
mass (see e.g. iLia et al.ll2002l ). In our simulation, the for- 
mation of stars is not accompanied by any kind of feedback 
processes. We have run tests with 50 times higher values 
of po.afr and/or 10 times lower values of Tfloor and did not 
find any major changes to our results. Both changes tend to 
produce dynamically hotter stellar disks as a lower cooling 
floor leads to stronger fragmentation of the disk and thus en- 
hances disk heating by dense substructures. A higher density 
threshold lowers star formation rates in the later, gas-poor 
evolution of our disk models thus diminishing the young, 
cold disk component. 



3 QUANTIFYING A DISK MODEL 

We start the analysis of our models by discussing in de- 
tail one model in halo A, which is named ARef. We choose 
this model as it produces a system which has certain struc- 
tural similarities to the Milky Way. The parameters applied 
are an initial radius of -Rgas = 60 kpc , a total mass of 
Mgas = 5 X lO^^Af©, a cooling time which is about twice as 
long as the free-fall time for the outermost gas (see figure 
[!}, a rotational velocity which increases from 120 kms~^ in 
the centre to 180 kms~^ at -Rgas and an orientation of the 
baryonic angular momentum vector parallel to that of the 
dark matter within 300 kpc. This orientation is close to the 
minor axis of the potential of the inner halo. The depen- 
dence of the final disk properties properties on these initial 
parameters is discussed in the following section. We split 
this section into subsections that deal with disk structure, 
angular momentum loss, kinematics and disk heating. 



3.1 Structural evolution 

Analyzing ARef at z — 0, 86% of the baryonic mass has 
formed stars and the remaining gas is situated mainly in 
a thin, extended gas disk (see top row in figure |8]). This 
;as fraction is in reasonable agreement with z '^ spirals 
Evoli et al.ll2011^ . As can be seen in the top left panel of 
figure [2] the star formation rate (SFR) rises steeply to a 
peak at 2: ~ 1 before falling by two orders of magnitude by 
z = 0. 

In the middle left panel of figure [2] we show the evolu- 
tion of the stellar surface density profile E(i?) of model ARef 
from 2 = 1.1 until 2 = 0. It can be seen that the first stars 
form in a concentrated bulge-like structure at i? < 2 kpc. 
This structure is actually triaxial because the infalling gas 
adjusts to the potential imposed by the halo. Its surface 
density only increases by a factor of ~ 2 by 2 = 0. The 
surface density profile outside the bulge region is already 
exponential in this early phase. Once most of the baryons 
have cooled, the infalling gas and the forming stars settle 
in a well-defined disk, which grows transforming the cen- 
tral potential into an axisymmetric configuration. The disk 
grows inside-out showing a t runcation at a radius -Rmax(2 ) 
Eis observed in disk galaxies (|van der Kruit fc Searle|[l98ll ). 
This is a consequence of high-angular-momentum gas liv- 
ing initially at greater radii, lower densities and thus longer 
cooling times. Moreover, according to our assumed star for- 
mation law, star formation timescales are longer at lower 
surface densities. The break in E(_R) moves outward with 



time to around i? ~ 19 kpc at 2 = 0. There is observa- 
tional evidence for inside-out formation of real disk galaxies 
(|Wang et al.ll201l|) and for o utward-movement of their trun- 
cations JAzzolini et aLlbOOSl ). 

Due to the numerically imposed cooling fioor Tfloor and 
the density threshold pthrosh, there is an imposed star for- 
mation threshold surface density Ethrosh ■ As is shown in the 
first row of figure |8] the gas disk at 2 = extends to larger 
radii, where this threshold is not reached. This also indicates 
that the truncation of the initial spherical distribution does 
not play a role here. The extende d gas d i sk is also slightly 
warped, as in real disk galaxies (jSancisil 119761 ). The warp 
is a result of a slight misalignment between the disk and 
the late-infall gas (see sections 5 and 7 for a discussion of 
misalignments) . 

Exponential disk profil es with central conce ntrations 
are common in real galaxies ijde Vaucouleurslll958l ) and also 
in simulations (jKatz fc Gunnlll99ll ), although most simu- 
lations significantly over- predict the mass in th e central 
bulge component (CS09, IScannapieco et al.l I2OI0I . but see 
iBrook et al.l I2OI2I for simulations of lower mass (nearly) 
bulge- less disks). The E(7?) profile at 2 = as displayed in 
the middle right panel of figure [2] is well-fitted by a profile of 



the form Eo,i exp ( — -n^ ) + So, 2 exp 



f_JL 



(*) 



(^) 



which 
The 



combines an exponential disk with a Sersic profile 
Sersic-indices found for our disk models are n ~ 0.5 — 2., sim- 
ilar to observed pseudo-bulges. The scale-lengths for ARef 
are Ri, = 1.0 kpc and Rd — 2.9 kpc. Calculating a disk-to- 
bulge mass ratio from this profile alone yields D/B ^ 12. 

In the bottom left panel of figure [5] we display the ver- 
tical stellar density profile p{z) for all stars in model ARef 
that are located in an annulus of width 1 kpc centred around 
R = 5 kpc at 2 = 0. The profile is well fitted by a double 
exponential with scale heights 420 and 1300 pc. The vertical 
profile does not vary significantly within R < 10 kpc, but 
at larger radii the disk flares and the scale height increases 
by a factor of ~ 2 until -Rmax- Thus the outer vertical struc- 
ture of this model is in conflict with observations of con- 
stant vertical structure withi n the exponential disk region 
(Ivan der Kruit fc Searlelll981^ . In the following section we 
will show that this conflict is absent for some of our models. 
The flaring and the 'thick' disk are caused by old stellar ma- 
terial formed in the bulge- formation period around 2 ~ 1. 
From then on the scale-heights flrst become smaller due to 
the growth of a young, cold disk before beginning to increase 
again due to disk heating. This increase is stronger in the 
inner disk as the outer disk has a younger mass-weighted 
age. The ARef disk thus has similar structural character- 
istics to the Milky Way disk (jjuric et all I2OO8I : iMcMillanI 
I2OIJ) , although its disk is ~ 40 % thicker. 



3.2 Angular momentum loss of infalling gas 

The radial structure of the disk discussed above is deter- 
mined by the potential and the angular momentum content 
of the stellar population. The concentrated peak in E(r) 
produced in the early stages of star formation indicates an 
unexpectedly high amount of low angular momentum mate- 
rial and thus angular momentum losses. We investigate this 
in figure [S] in which we plot the evolution of the total spe- 
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Figure 2. Properties of model ARef. Top Left: Star formation rate as a function of redshift z. Top Right: Stellar kinematics as a function 
of disk radius R at z = 0. Only stars within 500 pc of the midplane were taken into account. In purple we overplot the spherically- 



averaged circular velocity V^ir 



GM{<r) 



Middle Left: The evolution of the surface density profile S{-R) from z 



Middle Right: Surface density profile S(iJ) at z=0 fitted with a two-component profile Eo,l Gxp { — -ff- 



1.1 to z = 0. 
with 



.Rb) 



■ So, 2 exp 

scalelengths R^ = 2.9 kpc and Rf, = 1.0 kpc and n = 0.47. Bottom Left: Normalized z = vertical density profile at i? = 5 kpc fitted by 
a (double-)exponential with /ithin = 420 pc and /ithick = 1300 pc Bottom Right: Circularity distribution for all stars at z = 0. 



cific angular momentum of populations of baryons initially 
located in concentric shells of 5 kpc width. The material ini- 
tially located at _R < 10 kpc loses >~ 99% of its angular 
momentum within the first 500 Myr. The material in the 
next two shells still loses >~ 90%, the outermost material 
~ 50%. The strongest losses occur in the phase when the 



central potential is transformed from a triaxial to a disk-like 
configuration, in the process of which the bulge forms. More- 
over, the losses are subsequently delayed for outer shells, in- 
dicating that the angular momentum is lost as the gas cools 
and moves to the centre. 
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(Aunier & Binney 2009), whereas the thick disk shows ~ 5 
(Roddy ct al. 2006). To be able to compare different models 
better, we sort the stars according to their age, bin them in 
equal mass bins, and sort those according to the fraction of 
the z = stellar mass that has formed up to the time in con- 
sideration. We do this for five different redshifts. The black 
2 = line shows that Kot/f^ monotonically decreases with 
increasing age. Only the youngest stars are compatible with 
thin disk kinematics. The oldest ~ 1/5 is dispersion domi- 
nated and corresponds to the bulge peak seen in the bottom 
left panel of figure [2l Most of the galaxy thus has thick disk 
kinematics. Considering the evolution of Vrot/o'z over time, 
we can see that the bulge kinematics do not change from 
2 — 0.9 to 2 = 0.0, whereas the disk populations show de- 
creasing Kot/o"z ratios over time, i.e. the disk populations 
are heated, indicating that at least 1/3 of the stellar particles 
are born with thin-disk-like kinematics. 



Figure 3. Angular momentum loss as a function of time for model 
ARef. The initial gas sphere has been divided into 12 spherical 
shells of 5 kpc width each. For each shell the absolute of the total 
angular momentum vector j; of all particles initially in shell i is 
calculated at each output. 



3.3 Kinematical properties 

We move on to analyze the kinematics of ARef. The top right 
panel of figure [2] depicts the rotation velocity VJ-ot and the 
velocity dispersions an, a^ and a^ of the disk stars and the 
spherically averaged circular velocity Vcirc(7') = -v/GA/(r)/r 
profile at 2 = 0. Kirc peaks at ~ 250 kms~^ at i? -^ 10 kpc 
and gently decreases outside. The stellar rotation curve 
shows a broad peak at 7 < _R/kpc < 15 at Kot ~ 240 kms~^. 
For velocity dispersions, we find a decreasing profile out to 
7? ~ 12 kpc and an > o^ > cjz in the rotation-dominated 
region. The kinematical properties of the ARef disk are thus 
also similar to , but slightly hotter than the Milky Way disk 
JAumer fc Binnevll2009l : ISchonrichll2oT^ '). 

As a tool to quantify disk-to-bulge ratios, the circularity 
e = Vrot lViot.mB.^{ E) for a stellar pa rticle at energy E has 
been widely used (JAbadi et al.ll2003r ). Particles on a perfect 
circular orbit have e = 1 and thus a disk will show a peak 
close to e = 1, whereas a non-rotating bulge will have a 
peak at £ = 0. The bottom right panel of figure [2] depicts 
the e-histogram for ARef, which shows a distinct peak at 
e = 0.95 and only a shallow peak at e = 0, corresponding to 
a D/ B ~ 4. Thus the kinematic estimate for D/ B is signifi- 
cantly lower than the one from S(_R)-fitting , consistent with 
e.g. the results of IScannapieco et al.l l|2010f) for simulations 
of disk formation from cosmological initial conditions. 

However, this analysis is not capable of distinguishing 
between old and young stars and is also not good at quan- 
tifying how thin and dynamically cold a disk is. In the left 
panel of figure [4] we therefore employ a plot of the rotation- 
to-dispersion ratio Kot/o"z as a function of the age of the 
stars. We define a disk plane perpendicular to the total stel- 
lar angular momentum at each output time t. Kot is cal- 
culated as the mean tangential velocity and a^ as the rms 
velocity in the direction of the rotation axis of all stars in 
a certain age bin. Keep in mind that the Milky Way young 
thin disk has Kot/o"z >~ 20, the older thin disk has ~ 10 



3.4 Disk heating 

In order to analyze disk heating as it occurs in ARef in 
more detail, we plot the evolution of az for several coeval 
populations through time in the right panel of figure (4] 
The stars that form earliest and most centrally (see mid- 
dle left panel in figure [2| show high, almost isotropic veloc- 
ity dispersion (top right panel in figure [2]). These stars are 
strongly heated within a few 10* years (black line) or born 
hot (red line) . Within ~ 1 Gyr after the first star formation 
in the models the initial Gz has already become significantly 
lower, indicating the transition from bulge to disk forma- 
tion. Within another ~ 0.5 Gyr az at birth has dropped 
below 20 kms~^. Only the stars that form latest form with 
a smaller (Jz ~ 10 kms~^. As will be shown below, the high 
initial velocity dispersions of the first stars in this model are 
an effect of the structure of the dark halo. Test simulations 
in spherically symmetric haloes do not show this effect, but 
a significantly reduced initial velocity dispersion of the stars 
that form first, of the order ~ 10 — 20kms~^, as seen for the 
young populations in model ARef. The velocity dispersion 
of stellar particles in numerical simulations is limited by the 
applied star formation thresho ld, resolution and modeling of 
gas physics (JHouse et al.ll201ll ) and the formation of a real- 
istically cold young thin disk population is thus not possible 
in our current models. 

All the disky coeval populations show an increase of 
(Jz over time. Such continuous disk heating has been in - 
ferred for solar neighbourhood stars (JHolmberg et al.ll2009r ) 
and is usually linked to scattering of stars off disk substruc- 
ture such as transient sp irals and/or giant molecular clouds 
(jJenkins fc Binnevlll990l ). For our models we cannot, how- 
ever, exc lude that this effect i s partly due to numerical heat- 
ing (e.g. ISteinmetz fc Whitelll997l ). A closer analysis of the 
phenomenon over all of our models reveals that the strongest 
heating always occurs when distinct substructure in the form 
of a bar or spirals is present. A look at the evolution of Oz 
at different disk radii R over time reveals that due to the 
combined effect of the heating of coeval populations and the 
inside-out formation of our models in the inner disk regions, 
az first drops due to the formation of a thin disk and then 
increases due to disk heating, whereas at outer radii the for- 
mation phase lasts longer and az remains almost constant 
after an initial drop. 
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Figure 4. Disk Heating in model ARef: Left: Tiie ratio of tiie rotation velocity of stars to their vertical velocity dispersion Vrot/o"z as a 
function of the fraction of stellar mass that has formed until the formation time tfavu\ oi a star particle relative to the total stellar mass 
at z=0 for model ARef. The first stars are found at x-value 0, the ones formed at 2: = at x-value 1. Curves are plotted for redshifts 
z = 0.9, 0.6, 0.4, 0.2 and 0. The endpoints of curves indicate the fraction of stellar mass that has formed until the given redshift, e.g. 
82 % for 2 = 0.6. The evolution of Viot/o'z over time can be inferred by considering the various colored lines at one specific x-value. At 
each output redshift z we define a disk plane perpendicular to the total stellar angular momentum. Vi-ot is then calculated as the mean 
tangential velocity and ctz as the rms velocity in the direction of the rotation axis of all stars in a certain age bin. Right: Vertical velocity 
dispersion Cz (t) as a function of time for coeval populations in model ARef. The stellar population is binned into 50 equally spaced bins 
in redshift z, a representative selection of which is depicted by the various colors in the figure. The starting point of a line indicates the 
formation time of the population in consideration. 



3.5 Summary 

In conclusion, our setup is capable of producing disk galaxy 
models with realistic structural and kinematic properties. 
Disk formation proceeds inside-out and follows an early 
bulge-formation episode, during which the central halo po- 
tential is transformed into a disk-like configuration. This 
transformation is connected to angular momentum losses 
for the infalling gas. Moreover, a thick disk component in 
the vertical profiles originates from this epoch. We confirm 
that for a disk/bulge decomposition, a kinematic analysis is 
superior to a surface-density decomposition. However, cir- 
cularity distributions alone, as often applied, are not ideally 
suited to distinguish properly between thin and thick disk 
components, which is why Vrot/crz ratios should be used in 
addition. Our model does show a significant thin disk com- 
ponent, however its Kot/fz is limited to values ~ 10 com- 
pared to 20 for real young thin disks, due to the numerical 
lower limit we impose on CTz . The disk undergoes secular disk 
heating and shows a significant decline of zero-age velocity 
dispersion of stars. 



10.0r 




1.2 1.0 0.8 0.6 0.4 0.2 0.0 



Figure 5. Star formation rates as a function of redshift z for 
various model. 



4 HOW TO MAKE BETTER DISKS 

In this and the next section we discuss the dependencies of 
our models on the parameters that define them: 

(i) The mass of the model Mgaa, the initial radius -Rgas 
and thus the initial density profile p(-R) 

(ii) The rapidity of cooling and thus accretion of gas onto 
the forming galaxy represented by the cooling time normal- 
ization tcool.O 



(iii) The initial angular momentum distribution of the gas 
represented by the rotation velocity profile Kot(^) 

(iv) The orientation of the angular momentum vector of 
the gas (see also the following section) 

Our goal is to understand which physical conditions favor 
thin disk formation in ACDM haloes. An overview over our 
models and their parameters can be found in Table[T] In this 
section, we will mainly focus on halo A. For the parameters 



Mrrs^s, Rs 



tc 



1,0 and Vrot(R) we could not find significant 
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Figure 6. Angular momentum loss for different models. The left panel features different orientations with otherwise constant parameters, 
whereas the right panel focuses on different representations of physics at the same orientation. The definition of angular momentum loss 
is as in figure [3] The solid lines represent all baryonic particles in the simulations, the dashed lines only the particles that initially reside 
at _R < 0.2 flgas. 



halo-to-halo differences. Most differences because of orienta- 
tions are discussed in Section 5. As in the previous section, 
we split our analysis into subsections that deal with star for- 
mation, angular momentum loss, disk structure, kinematics 
and disk heating. 



4.1 Gas infall and star formation 

We begin with a discussion of how the model parameters 
influence the formation histories of our models. The star 
formation rates of the models displayed in figure [S] reveal 
that the process of gas inflow depends mainly on tcooi.o and 
the loss of angular momentum, which is put into the initial 
models via Kot(-R). Models ARef (black line), ARef-|-Maj 
(purple) and AMaj2C (turquoise) share the initial density 
profile and cooling time. All models show a steep rise fol- 
lowed by a softer decrease in SFR. ARef-|-Maj and ARef 
differ only in orientation, consequently they show qualita- 
tively similar SFRs. ARef-|-Maj loses a smaller amount of 
angular momentum (see below), thus forms a more extended 
structure with lower densities and has longer star forma- 
tion timescales, a lower peak SFR and a higher late-time 
SFR. The second peak for ARef-|-Maj corresponds to star 
formation in a ring structure (see the following section). 
ARef-l-Maj and AMaj2C differ only in Kot(i?) with AMaj2C 
having significantly less angular momentum, which leads to 
a much more concentrated object, higher densities, higher 
cooling rates, higher SFRs, an earlier peak and a steeper de- 
crease in SFR. The same is true for a comparison of AExtlC 
(red line), AExt (green) and AExtLoAM (blue). They are 
all more extended with a higher initial gas mass than the 
models mentioned so far and thus represent a different ini- 
tial density profile. AExtlC cools 4 times faster than AExt 
and AExtLoAM has 3 times less angular momentum than 
AExt. In both cases this results in a higher and earlier peak 
in SFR. 

Directly connected to the SFR is the remaining gas frac- 



tion in the disks. At z = these are almost independent of 
cooling time, as the cooling and star formation timescales 
are small compared to the simulation time and the gas disk 
outside i?max contains a negligible amount of mass. Thus 
AExtlC and AExt both show a final gas fraction of 18%. An- 
gular momentum variation can significantly alter the mass in 
the outer gas disk with very long cooling and SF timescales. 
Thus, ARef+Maj (36%) and AExt (18%) have much higher 
gas fractions than AMaj2C (3%) and AExtLoAM (4%). Re- 
alistic gas ratios are 5 to 25% (lEvoli et al.ll201ll ). so that 
only our most extreme models are inconsistent with obser- 
vations. 



4.2 Why does infalling gas lose angular 
momentum? 

As we have shown in Section 3.2, gas in our models loses 
angular momentum as it falls to the centre. Figure [6] al- 
lows a closer inspection of angular momentum loss in sev- 
eral models. For each model the solid line represents the 
absolute value of the total specific angular momentum of 
all baryonic particles, whereas the dashed line represents 
only the particles that initially live within R < 0.2 Rgas- 
The left panel focuses on models which differ only in ori- 
entation and apart from AMaj2C, share the same density 
and velocity profiles and cooling times. Apart from ARef- 
Med and AMaj2C the central parts show higher angular 
momentum losses. The difference between central and total 
mass is biggest for model ARef as discussed above. Focus- 
ing on total angular momentum losses, the models aligned 
with the major halo axis (ARef-|-Maj and AMaj2C) show 
the lowest losses. They are about 10% at t = 8 Gyr and 
subsequently only increase for ARef-|-Maj which develops 
strong misalignments between inner and outer galaxy, lead- 
ing to higher losses than in AMaj2C, which forms a well- 
aligned disk. ARef45, which starts with its orientation 45 
degrees offset, has lost 40% of its total angular momentum 



© 2012 RAS, MNRAS 000, [T]-?? 



10 M. Aumer and S.D.M. White 



150- 



100 



ARef 

ARef180 

ARef+Med 

ARef-Med 

AExt180 



50- 




t [Gyr] 

Figure 7. Mass-weighted reorientation angle of all baryons as 
a function of time for models ARef-Med, ARef-|-Med, ARef, 
AReflSO, AExtlSO. The angle is measured in a fixed coordinate 
system between the initial angular momentum vector orientation 
at 2 = 1.3 and its given orientation at time t. 



at t — 8 Gyr. The rest of the models, which all started ori- 
ented perpendicularly to the minor axis, show losses well in 
excess of 50%. Only a small part of these l osses, perhaps of 
the or der of 10% can be of numerical origin. iKaufmann et al.l 
l|2007h reported 10% losses for their highest resolution model 
in a spherical halo, the resolution of which is similar to our 
models. 

As halo A has an approximately prolate shape 
l|Vera-Ciro et al.ll2011f ). the axis ratio of the halo potential 
in the plane perpendicular to the major axis is close to 1, 
i.e. almost axisymmetric, whereas it is ~ 0.7 if the disk axis 
is perpendicular to the major axis. Thus the formation of 
a disk in the former case requires a less severe transition 
of the central halo potential into an oblate, axisymmetric, 
disk potential, than is the case for the latter models. The 
strongest angular momentum losses in our models occur in 
these initial transformation phases. As our spherical initial 
gas distributions are not in equilibrium with the gravita- 
tional potential of the halo, the gas, which is significantly 
less dense than the dark matter in the region in consider- 
ation, adjusts to the triaxial potential in the early phases 
of our simulations. Due to its rotation the shape of the gas 
distribution becomes offset from the shape of the potential 
if the halo symmetry axis and the axis of rotation do not 
agree. Gravitational torques then lead to the transfer of an- 
gular momentum to the halo. This explains the initial loss 
of angular momentum in all models expect those oriented 
parallel to the major axis, as seen in figure |6] The differ- 
ences in the evolution of models at the same orientation or 
with orientations perpendicular to major, which all show 
similar axis ratios of the equipotential ellipses in the plane 
defined by their initial orientation, appear at later times (af- 
ter ~ 0.5 — 1 Gyrs) and coincide with phases of reorientation 
of the angular momentum vector (see below) . 

In figure [T] we illustrate this reorientation of the total 
baryonic angular momentum vector in models ARef-Med, 
ARef-|-Med and ARef. Strong reorientation starts to develop 



after ^ 0.5—1 Gyrs. At these times, the assembly of baryons 
in the centre of the halo has led to a transformation of the 
central halo potential shape, which apparently can induce 
strong reorientation of the baryonic angular momentum vec- 
tor. ARef shows the lowest amount of reorientation and has 
the lowest losses of angular momentum. ARef-Med shows 
continuous reorientation, whereas ARef-|-Med settles into a 
new orientation relatively quickly. Consequently ARef-Med 
loses the most angular momentum (see figure [S}. 

In the right panel of figure [6] we focus on models with 
the same rotation axis. AReflSO, which differs in orientation 
by 180 degrees from ARef, and loses more angular momen- 
tum, as it shows strong reorientation, (figure [7]). However, 
the losses are smaller than in ARef-Med, showing that an- 
gular momentum losses depend not only on reorientation 
and on the shape of the potential in the disk plane, but 
also on the complex details of the interaction of hot gas, 
dark halo and stellar component. AExtlC, which is more 
extended, more massive and has a higher total angular mo- 
mentum content, shows a smaller total angular momentum 
loss. AExt, which differs from AExtlC only by a 4 times 
longer cooling time, loses approximately the same amount 
of angular momentum, but on a longer timescale, illustrat- 
ing again that the loss occurs as the gas cools to the centre. 
AExtLoAM, which has the lowest amount of angular mo- 
mentum and goes bar unstable (see figure [SJ, shows smaller 
losses than AExt, which only differs in Vrot{R)- 

For fixed values of other parameters, shorter cooling 
timescales, stronger reorientation and higher initial densi- 
ties lead to increased angular momentum loss. Higher initial 
angular momentum also leads to greater fractional losses. 

4.3 Structural properties 

We continue with a study of the infiuence of the model pa- 
rameters on the structure of the disk models. In figure]^ we 
analyze S(r) profiles for several models in halo A. Here we 
focus on the models, that produce the strongest disks. The 
others are discussed in the following section. Model AExtlC, 
which has more initial angular momentum, a higher mass 
and a larger initial radius than ARef (see bottom left panel 
in figure [2J produces a more extended disk with a disk-scale 
length Rd = 5.3 kpc and no truncation out to 30 kpc. The 
four times longer cooling time of AExt mildly extends the 
scale-length to 6.1 kpc, but produces a truncation around 
20 kpc as the outermost material cools too slowly to form 
stars. The bulge profiles are very similar in the two models, 
showing that the formation of a centrally concentrated com- 
ponent due to transformation of the central halo potential 
depends little on cooling time. 

Model AExtLoAM goes bar unstable and thus shows a 
much stronger central component with a longer scale-length 
(1.2 vs 0.9 kpc). The bulge profile extends to 8 kpc. The pro- 
file outside can be approximated by an exponential disk, but 
clearly shows a ring structure, that is also visible in figure 
[S] in both stars and gas. Pseudo-bulge and ri ng are results 
of bar -induced gas flows as discussed e.g. in lAthanassoulal 
(■19921). 

In model AMaj2C there is no severe initial central trans- 
formation. The profile can be fitted by a single exponential, 
that, because of the very low angular momentum losses, 
has a larger scale-length than ARef (4.5 vs 2.9 kpc) de- 
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Figure 8. Edge-on and face-on surface density projections for stars (left panels) and gas (right panels) at 2: = for models (top to 
bottom): ARef, CExtCosmo, AExtLoAM, ARef-Min, CExt-hAM and C-Maj. 
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Figure 9. Radial stellar surface density profiles, S(iJ) at z = 
for models CExt+AM, AExtlC, AExt, CExtCosmo, AExtLoAM, 
C-Maj, AMaj2C and ARef-Min. Overplotted are fits of the type 
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Figure 10. Vertical density profiles of models CExt+AM, C- 
Maj, ARef-Min, AExtLoAM, AExt and CExtCosmo. If sensible, 
exponential fits are overplotted and the scaleheights are noted. 



spite the low initial angular momentum content. The low 
angular momentum also leads to a small truncation radius 

-Rmax ~ 12 kpc. 

The coldest disks form in AExt and CExtCosmo, which 
share initial density and angular momentum profiles as well 
as cooling timescales, and which both undergo little reorien- 
tation. They show small halo-to-halo differences in terms of 
T,{R). AExt has a more massive bulge component, a slightly 
less extended disk (R^ = 6.1 versus 7.0 kpc) and a slightly 
smaller truncation radius. This underlines that the effects of 
the parameters discussed in this section are very similar in 
the two cases. 

In figure[Tn]we compare vertical profiles of several mod- 
els at _R = 5 kpc. In contrast to the double-exponential of 
ARef (see figure [2]), AExt shows a single exponential profile 



with a slightly reduced exponential scale-height h^ — 400 pc 
(compared to 420 pc for ARef). We also present the profile of 
CExtCosmo. It has an even thinner disk with hz — 330 pc, 
consistent with the reduced bulge fraction discussed above. 
Moreover, the flaring of these disks is significantly reduced 
compared to ARef. Scale-heights for CExtCosmo are h^ < 
400 pc at i? < 10 kpc and h^ < 550 pc out to i? ~ 20 kpc. 

The edge-on and face-on surface density projections in 
gas and stars of figure [8] (second row) nicely reveal that 
CExtCosmo sustains spiral substructure in the stellar and 
gas disk until z = 0, unlike ARef, which at this time has 
a significantly lower SFR. A comparison to ARef also high- 
lights the less prominent bulge structure and the lack of an 
underlying thick component. A longer cooling time and a 
higher angular momentum content result in a prolonged star 
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Figure 11. Circularity distributions for models ARef, ARef-Min, 
AExt, AExtLoAM, AReflSO and C-Maj. 



formation in the disk of CExtCosmo. Together with a small 
amount of reorientation of the angular momentum vector, 
this leads to the continuous existence of a relatively thin 
stellar disk preventing the formation of a thick component. 
In comparison to these extended disks the low angu- 
lar momentum model AExtLoAM, which goes bar unsta- 
ble shows a thick vertical profile at i? = 5 kpc (in the bar 
region), which is fit neither by an exponential nor by an 
isothermal profile. In the edge-on view in the third row of 
figure |8] the peanut-shape of the stellar system is clearly 
visible. 



4.4 Kinematical Properties 

The significant parameter dependencies found above are also 
reflected in the stellar kinematics of our models. In terms 
of the circularity distribution, longer cooling times and the 
resulting thinner and more extended disks result in a slightly 
enhanced disk peak and a reduced bulge peak as may be 
seen by comparing AExt and ARef in figure [TT] The massive 
bar in AExtLoAM significantly diminishes the e ~ 1 peak 
and shifts material into a wide, asymmetric distribution at 
-l<e< 0.8. 

The parameter dependence of the disks is summarized 
in the rotation-to-dispersion-ratio plots of figure [121 In the 
upper panel we compare models ARef, AExtLoMa and 
AExt2C. Compared to ARef, AExtLoMa is more extended 
(^gas ~ 100 vs 60 kpc). This leads to an increase of Kot/fz 
for all generations of stars, a reduction of the bulge- fraction, 
but no significant increase in the thin disk population. A 
lower density, which already results in slower cooling (cf. 
figure [T|, thus helps reduce the impact of initial bulge forma- 
tion. An increase in mass by 50% has little impact as may be 
seen by comparing AExt2C to AExtLoMa. It leads to more 
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Figure 12. Rotation-to-dispersion ratio as in figure |4l but now 
for different models at z = 0. Top: The effect of the initial ra- 
dius ijgas and mass Mgas illustrated by models ARef, AExtLoMa, 
AExt2C and AExtHiMa. Middle: The effect of the normalization 
of the cooling time tcool,o illustrated by models AExtlC, AExt2C 
and AExt. Bottom: The effect of the angular momentum content 
/ rotation velocity profile Vrot{r) illustrated by models AExt, 
AExtMedAM and AExtLoAM. 
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Figure 13. Vortical velocity dispersions ^^(i) as a function of time for coeval populations in models AExtLoAM (left) and CExtCosmo 
(right). See figure |4] for an explanation of how these plots are constructed. 



z — thin disk stars as it increases the surface density and 
thus prolongs star formation in the disk outskirts. Increasing 
the mass again by a factor of 2, as in AExtHiMa, reduces 
the rotation-to-dispersion-ratio significantly. The bulge frac- 
tion is reduced due to the longer cooling time applied (see 
middle panel). The disk goes bar unstable, however, and the 
bar subsequently heats the disk. 

The second panel compares the models AExtlC, 
AExt2C and AExt and thus reveals the impact of the cool- 
ing time. Clearly a longer cooling time leads to an increase 
of Viot/cTz for all generations of stars. The increase is espe- 
cially significant for the old population. AExt only shows 
~ 5% bulge stars. From figure [4] we know that Kot/o"z de- 
creases due to disk heating as stellar populations age. As 
we know from figure [5] that AExt has a much younger pop- 
ulation than AExtlC, the difference for the disk stars can 
be attributed mainly to the younger ages of the population. 
Kot/o"z for the thin disk population is similar for all models, 
confirming that this is mainly determined by resolution and 
by the prescriptions for star formation and for cooling. 

The bottom panel of figure[l2]concerns the initial angu- 
lar momentum content of the models. Compared to AExt, 
Kot(-R) in AExtMedAM has been reduced by 1/3 and by 
2/3 in AExtLoAM. Interestingly, AExtMedAM shows a 
slightly smaller dispersion-dominated component. This indi- 
cates that the transfer of angular momentum from baryons 
to the dark halo plays a crucial role in the initial phase of 
bulge formation/transformation of the central potential. For 
the younger population, AExt is colder than AExtMedAM, 
as it is more extended and thus has longer star formation 
timescales for the outer populations. AExtLoAM has been 
discussed above and shown to go massively bar unstable. 
The consequence is a drastic reduction of Vrot /ctz for all 
ages. 



4.5 Disk heating 

To illustrate differences in disk heating, we depict in figure 
1131 the evolution of a^ (t) for several coeval populations of 



stars in models AExtLoAM and CExtCosmo. Compared to 
model ARef, presented in figure [J] they both have lower ini- 
tial densities and slower cooling times. Like ARef they show 
continuous disk heating for all populations. Moreover, the 
oldest populations undergo strong heating in the first several 
10* years, but, unlike ARef, the slower initial transforma- 
tion of the central potential prevents hot initial dispersions 
o"z,ini > 50 kms~^. 

The low angular momentum model AExtLoAM shows 
the strongest heating of old components (to a^ ~ 
120 kms~^) due to the most concentrated surface density 
profile. At i ~ 11 Gyr after the massive bar has developed 
(cf. figure [8)l a strong enhancement in the heating of disk 
populations is visible. In less than a Gyr, az{t) increases by 
about 50% for all populations except the oldest, which are 
dispersion-dominated and are hardly affected by the bar. An 
inspection of radial velocity dispersion aR{t) reveals that in- 
plane heating is already enhanced at t ~ 9 Gyr, when the 
bar starts to form. The bar also increases the birth veloc- 
ity dispersions of stars by about 50%. The final bar/disk 
populations are heated more efficiently and are thus signifi- 
cantly hotter than in ARef. Analyzing az{t) for radial bins 
shows that the increase occurs at all radii less than ~ 8 kpc 
and also in the central kpc. The stellar populations outside 
R >~ 10 kpc (outside the bar region) are however not af- 
fected. 

For the disk populations in the high angular momentum 
model CExtCosmo, az{t) behaves as for the correspond- 
ing populations in ARef, continuously increasing with the 
birth dispersions decreasing with time. The initial heating 
for the oldest populations is however significantly reduced 
compared to ARef. Consequently the transition between the 
oldest and the subsequent populations is not as strong as 
in the other models. This is reflected by the increase in 
D/B for disk models with longer assembly timescales. Still, 
CExtCosmo and all other disk models show a decrease in the 
birth velocity dispersions with time. This behaviour is a con- 
sequence of the triaxial halo and its substructure, and does 
not occur in simulations within idealized spherical haloes. 
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4.6 Summary 

In this section, we have studied the influence of coohng 
timescales, angular momentum content, initial density pro- 
file and mass on our models. We show that the formation 
timescale increases with increasing cooling time and increas- 
ing angular momentum content. The angular momentum 
loss in our models increases mainly with increasing initial 
ellipticity of the dark matter potential in the disk plane. If 
reorientation of the baryonic angular momentum vector oc- 
curs in a model, it adds to the loss of angular momentum. A 
slower formation process and thus a slower transformation 
of the central potential yields thinner and more dominant fi- 
nal disks and weakens the formation of prominent bulge and 
thick disk components. Massive bulges can only be avoided 
if the initial potential contours in the disk plane are almost 
circular. Double-exponential vertical disk profiles are sup- 
pressed by a slow formation timescale. All our disks show 
continuous secular heat ing due to (spiral) substructure (e.g. 
I Jenkins fc Binnevll 199(1 ) . so that younger disks are thinner. 
Disk heating is st rong l y enh anced by bars, as was recently 
also discussed by, S aha e t al. (2010'). Additionally, heating is 
strongest for the old bulge populations and the birth velocity 
dispersion of stars in our simulations always decreases with 
tim e, which adds to th e age- velocity dispersion relation (see 
also lHouse et al.ll201ll ). Higher angular momentum content 
leads to more extended disks and larger truncation radii. 
The latter also increase with faster cooling. Bars form, if, 
at fixed other parameters, the mass exceeds or the angular 
momentum falls below (due to loss or initial lack) a sta- 
bility limit consist ent with standard stability criterions for 
bar formation (see lBinnev fc Tremaine 2008 §6.3). (Pseudo- 
)Bulges are enhanced by bar-induced gas inflows. 



5 (RE) ORIENTATIONS OF DISKS 

So far we have paid relatively little attention to the influence 
of the initial orientation of the angular momentum vector 
of the rotating gas sphere. We have shown that it plays a 
crucial role in angular momentum loss of the infalling gas 
(figure [6]). Moreover, figure [7] revealed that this orientation 
is not fixed, rather different models can undergo very differ- 
ent amounts of reorientation. Model ARef and many others, 
which we used to study the influence of parameters, were ini- 
tially oriented parallel to the angular momentum vector of 
the dark matter within 300 kpc. This orientation was chosen 
as it yielded the coldest disks in several test simulations and 
showed an amount of reorientation that is similar to that of 
the dark halo itself from z = 1.3 to z = 0.0. In terms of the 
principal axes of the triaxial potential ellipsoid, the orienta- 
tion of ARef is close to the minor axis, but produces slightly 
better results than an exact alignment with this axis. 

In halo C, this setup failed. The dark matter angular 
momentum is poorly aligned with the principal axes and 
thus the plane defined by this orientation contains strong 
variations of the vertical gravitational force, preventing qui- 
escent disk formation. However, orientations close to the mi- 
nor axis proved to work well, underlining that the influence 
of the shape of the halo potential is stronger than the influ- 
ence of its angular momentum. This flnding is consistent 
with previous findings that disks in cosmological simula- 
tions preferentially align with the minor axis of the halo 



teail in et al.l |2005| ). In this section, we therefore explore 
which orientations in haloes A and C are capable of pro- 
ducing thin disks. We do so in subsections, first discussing 
(re)orientations in minor and medium orientation models 
and the effects of reorientation on kinematics. Then we dis- 
cuss major/intermediate orientation models, the structure of 
peculiar models and how reorientation shapes the potential. 

5.1 Models with orientations perpendicular to the 
major axis 

As both haloes are close to prolate at z — 1.3 we start 
by exploring orientations, which lie in or close to the plane 
perpendicular to the major axis of the halo potential. We 
illustrate the findings for halo A in figure [141 where in the 
left half we show the initial orientations of the angular mo- 
mentum of the gas sphere and in the right half we show the 
orientations of the angular momentum of all the stars that 
have formed by 2 = 0. The coordinate system is fixed and 
identical in both halves. Zero degrees is defined by the ini- 
tial orientation of ARef and the numbers in brackets behind 
the model numbers indicate the angles by which the orienta- 
tions are offset from the plane in degrees. Initially all these 
models are within 12 degrees of this plane. We show mod- 
els aligned with the dark angular momentum within 300 
kpc (ARef, AExt), 180 degrees offset from this (AReflSO, 
AExtl80), ahgned with the minor axis (ARef-fMin, ARef- 
Min) and with the medium axis (ARef+Med, AExt+Med, 
ARef-Med, AExt-Med) and aligned with the final orienta- 
tion of model ARef (ARefEnd). 

The final orientations of the galaxies in halo A show a 
clear crowding of models around the minor axis. All models 
with orientations between ARef+Min and ARef-|-Med end 
up within ~ 20 degrees from the reference model ARef and 
close to the plane perpendicular to the initial potential ma- 
jor axis. Figure [7] shows a continuous reorientation of a total 
of 30 degrees for ARef. ARef-|-Med however first turns by 
~ 90 degrees before also showing a small, continuous reori- 
entation. Thus these models all seem to adjust to the 'cor- 
rect' orientation and then evolve rather quiescently. Models 
with different parameters sharing the same initial orienta- 
tion show very similar final orientations as illustrated by 
models ARef/ AExt and ARef-fMed/ AExt+Med. 

All the other models however end up with significantly 
different orientations, and there is no clear correlation be- 
tween their initial and final orientation. ARefl80 seems to 
end up close to its original orientation, however it flips by 
almost 120 degrees before slowly returning. If dark matter 
influenced the evolution only through its global potential, 
there should be no difference between two orientations dif- 
fering by 180 degrees, which is clearly not the case, so that 
halo rotation and transfer of angular momentum between 
the baryons and the dark matter must play a role in the de- 
scribed processes. Comparing AExtl80 and ARefl80, which 
share initial orientation, we detect a final offset of more than 
90 degrees. A smaller, but significant difference exists also 
between AExt-Med and ARef-Med. In figure [3 one can see 
that the initial turn of ARefl80 by ~ 120 degrees is also seen 
in AExtl80, where this process takes longer. The model then 
stabilizes in this orientation. Model ARef-Med, which starts 
from a medium orientation, shows continual and strong re- 
orientation. 
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Figure 14. Visualization of reorientation in halo A for models with initial orientations (approximately) in the plane perpendicular to 
the major axis of the halo potential. Left are initial orientations at 2 = 1.3, right are mass-weighted mean orientations of the stellar 
systems at z = 0. The coordinate system is static and was fixed to be perpendicular to the major axis of the halo potential at z = 1.3 
and to have the initial orientation of model ARef at degrees. The Greek letters and colors code the models and the numbers in brackets 
indicate the angle between the orientation and the plane in degrees. AQg stands for the orientation of the angular momentum of the gas 
at ij = 30 — 100 kpc at 2 = in the simulation of CS09. 
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Figure 15. As figure [Ml but for models in halo C. The initial orientation of model CExt+Min is at degrees. AQ stands for the disk : 
halo C of CS09. 
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Figure 16. Rotation-to dispersion ratios as defined in figure|4]but 
now for models with difl'erent orientation. A negative ratio indi- 
cates a component that is counter-rotating with respect to the to- 
tal stellar angular momentum. Top panel: Models perpendicular 
to the major axis ending up in 'disk orientation' and disk-like ro- 
tation: ARef-l-Med, ARef, ARef-|-Min, ARefEnd and AExt. Mid- 
dle panel: Models perpendicular to the major axis not showing 
a continuous disk formation, but all showing significant reorien- 
tation and developing counter-rotating components: ARef-Med, 
AReflSO, ARef-Min, AExt-Med, AExtlSO. Bottom panel: Mod- 
els parallel to the major axis: ARef-|-Maj, ARef-Maj, AMaj2C, 
AMaj, AMaj 180 and model ARef45. 



In figure [TSl we illustrate the situation in halo C. Zero 
degrees is defined by the initial orientation of CExt-|-Min, 
which was aligned with the halo minor axis. CExt-Min is 
180 degrees offset, CExt-Med and CExt-|-Med are aligned 
with the medium axis and CExtEnd is the final orientation 
of CExt-l-Min. CExt-l-Arb and CExt-Arb resulted from an 
error in the setup but are useful for this analysis, AQ marks 
the orientation of the disk in the cosmological galaxy for- 
mation simulations of CS09. Finally, CExtCosmo uses the 
z = 1.3 cosmological orientation as an initial condition. We 
see that the cosmological disk is well aligned with the mi- 
nor axis, around which many models again crowd at z = 0, 
in this case all models with orientations between those of 
CExt+Arb and CExt-|-Min in the initial conditions. CExt- 
Min and CExt-|-Med indicate that, for halo C, the direction 
180 degrees offset from the 'best disk' orientation is also a 
preferred orientation a.t z = 0. Among the models presented, 
only CExt-Arb ends up in a significantly different orienta- 
tion. 



5.2 How reorientation affects kinematics 

In the previous subsection we have shown that the orienta- 
tion of the angular momentum of infalling gas can undergo 
significant changes. In the previous section we have already 
shown that this leads to loss of angular momentum. Here we 
study the effects of these processes on the kinematics of the 
models. 

In figure[T6]we show rotation-to-dispersion ratios for the 
halo A models ending up close to the ARef final orientation 
(top panel) and for the other models of figure [TJ] (middle 
panel). In both panels, the models with longer cooling times 
show higher Kot/o"z for the disky populations as was dis- 
cussed in the previous section. There is, however, a clear 
dichotomy between the two panels. All models in the up- 
per panel show a monotonic increase in Vrot/cz from old to 
young stars, indicating that these are continuously forming 
disks. It also shows that ARef has the coldest disk popula- 
tion of the models sharing the same parameters, presumably 
because it also shows the smallest amount of reorientation. 
However, ARef+Min and ARef+Med, which start exactly 
in the plane perpendicular to the major axis show a smaller 
bulge population than ARef and ARefEnd, which start with 
a small misalignment to this plane. 

The models depicted in the middle panel do not end up 
close to the preferred disk orientation and all show counter- 
rotating young populations indicating the infall of gas mis- 
aligned with the rotation of the stellar object at late times 
and none of these objects show thin disk populations of stars. 

Circularity distributions for AReflSO and ARef-Min are 
shown in figure [TT] The compact, barred disk of ARef-Min 
has lost its e '^ 1 peak and shows a wide, asymmetric dis- 
tribution peaking at e ~ 0.6. It has a strongly enhanced 
counter-rotating population compared to ARef. ARefl80, 
in contrast, shows only a mildly diminished disk peak 
compared to ARef, but features a distinct, small counter- 
rotating disk peak at e ~ ( — 1). 

In figure [17] we analyze the radial profiles of vertical 
velocity dispersion az (R) for several of the models just dis- 
cussed. Unlike ARef, which shows a realistically declining 
(Tz{R), there are several models showing an outward increase 
in Gz(R) after a drop at inner radii, most strikingly AReflSO. 
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Figure 17. Vertical velocity dispersions a-z{R) as a function 
of disk radius R in the models ARef+Med, ARef, AReflSO, 
ARef+Min, ARefEnd. 



According to figure[71 AReflSO shows the strongest reorien- 
tation during its evolution and according to figure [6] also 
the strongest losses of angular momentum of all models, 
whereas ARef shows only weak reorientation. This trend is 
confirmed by AExtlSO and ARef+Med, which are intermedi- 
ate between ARef and AReflSO in terms of ctz (R) flaring and 
reorientation. ARef-Med, which shows the strongest angu- 
lar momentum losses and continuously strong reorientation 
shows very high a^ (R) at all radii and mild flaring. This also 
strengthens the conclusion that these two processes lead to 
high velocity dispersions in the final systems, especially in 
the outskirts. 

In the top panel of figure [18] we plot rotation-to- 
dispersion ratios for all models in halo C that end up with a 
thin disk population of stars (V/az > 5 and V/az decreas- 
ing monotonically with age). CExtEnd and CExtCosmo, the 
orientations motivated by the cosmological model and the 
final orientation of CExt+Min, produce the coldest disks. 
Similarly to AExt in halo A, these models are slightly offset 
from the minor orientation. CExt+Min, the minor orienta- 
tion run, still has a better disk than CExt-Min, which starts 
and ends at an angle of ISO degrees. Unlike in halo A, this 
orientation still produces a disk. The medium orientation 
runs CExt-Med and CExt+Med, which also start and end up 
at an angle of ISO degrees produce slightly worse disks due to 
the effects of reorientation. CExt+Arb, which starts from an 
orientation about 120 degrees offset from CExt+Min, also 
produces a disk in the preferred final direction. CExt-Arb 
is the only model not following the trend. The orientation 
changes strongly and the model results in a thick disk sur- 
rounded by misaligned components, each resulting from a 
different formation phase as indicated in the lower panel of 
figure [TSl 
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Figure 18. Rotation-to dispersion ratios as defined in figure |4l 
but now for models in halo C. Top panel: Models producing disks: 
CExt+Min, C+Maj, CExt-Min, CExt-Med, CExt+Med, CEx- 
tEnd, CExtCosmo, CExt+Arb. Bottom panel: Other models in 
halo C: CExt+AM, CExt-AM, C-Maj, CExt-Arb. 



5.3 Model with major and intermediate initial 
orientations 

Until now, we have focused only on models that start 
(nearly) perpendicular to the major axis. As already dis- 
cussed in the previous section, AMaj2C, which starts aligned 
with the major axis, also produces a relatively small, thin 
and exponential disk. As can be seen in the bottom panel 
of figure 1161 the same is true for AMaj, which employs 
slower cooling, but produces only a slightly colder disk. It 
was shown in figure [S] that the infall and star formation 
timescales are determined by angular momentum and cool- 
ing time. For these models with low angular momentum, 
star formation is shifted to earlier times, and consequently 
the stellar populations are subject to disk heating over a 
longer period. As a result the final difference in az is small. 
AMajlSO, which uses the same parameters as AMaj, but 
is initially ISO degrees offset, produces a disk with almost 
the same characteristics as AMaj. All these models reorient 
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away from the initial major axis, but by different amounts, 
~ 50 degrees for AMaj2C and AMaj, which end only ~ 5 
degrees offset, and ^ fOO degrees for AMaj f 80. They do 
not end up near the minor axis, but do not seem to be af- 
fected much by reorientation. Probably this is because their 
formation timescales are small and they are compact. Mis- 
alignments typically happen for material that comes in late 
at large radii and is thus hardly affected by the potential of 
the inner disk. 

As alluded to in the previous section, the plane defined 
by the major axis orientation in almost prolate systems has 
an almost axisymmetric potential and is thus advantageous 
for our models compared to the minor axis orientation. How- 
ever, models ARef-|-Maj and ARef-Maj, which start with sig- 
nificantly more angular momentum than AMaj /AMaj 180 do 
not produce disk-like objects. This is also true for ARef45, 
which starts from an unfavorable orientation half-way be- 
tween minor and major axes, as can be seen in the bottom 
panel of figure 1161 The reason for this failure is misalign- 
ments that develop between gas infall at early and late times, 
destroying an initial disk. 

In halo C, the situation is similar. C+Maj, which starts 
in major axis orientation with the same parameters as AMaj, 
also produces a disk. It also undergoes reorientation away 
from the major axis and ends up at an angle ~ 30 degrees 
from the final orientation of CExt-|-Min. However, C-Maj, 
which has an initial gas angular momentum orientation 180 
degrees apart from C-I-Maj and a more extended initial ra- 
dius and thus higher angular momentum, does not produce 
a nice disk, showing outer misalignments as in ARef-|-Maj. 
This model shows that misalignments are common for initial 
orientations parallel to the major halo axis, also for medium 
angular momentum content. 

The model CExt-(-AM is initially aligned with the an- 
gular momentum of the halo at i? < 300 kpc, which is at 
an angle of ~ 45 degrees to the major axis. This is why the 
model experiences similar problems to ARef45. The same is 
true for CExt-AM, which starts at an angle of 180 degrees 
to CExt+AM. 



5.4 The structure of peculiar models 

Although our models were designed to produce disks, not 
all of them end up as such, as we have shown above. Here 
we briefly analyze the structure of three of those models in 
more detail. 

ARef-Min is one of the models, which start oriented per- 
pendicular to the major axis and show continuous, strong re- 
orientation, which leads to strong angular momentum losses 
and high velocity dispersion (see middle panel in figure 
I16p . We have included surface density projections of ARef- 
Min in figure [8] (fourth row) . Due to the angular momen- 
tum losses the resulting disk is compact, barred and thick 
with an exponential scale height of h^ = 1.3 kpc (see fig- 
ure [10]). Its exponential-like radial profile is very steep with 
Rd ~ 1.8 kpc as shown in figure [9] The galaxy has a high 
bulge fraction, D/B ^ 1. Even at 2; = the remaining gas is 
concentrated in the central bulge and does not show a disk 
structure. Also the stellar disk structure is mostly bar-like 
and is surrounded by a spheroidal star distribution. Its circu- 
larity distribution ffigurefTTj) shows a wide, asymmetric peak 
at e ~ 0.6. The system might be classified as a barred, fast 



rotating early-type galaxy applying a classification scheme 
as used for the ATLAS^-" survey (E msellem et al. 201l| ). 

CExt-|-AM represents the models starting with inter- 
mediate initial orientations. In the fifth row of figure [51 we 
show its 2 = surface density projections. The gas lives in 
a perfect disk, whereas the stellar disk shows two regions: 
a dominant bulge-like central structure, which is elongated 
perpendicular to the gas disk, and outer regions dominated 
by rings and separated from the centre by a drop in sur- 
face density around _R ~ 9 kpc (see figure [9]). The vertical 
profile at _R = 5 kpc can be fit by a thick exponential with 
hz ~ 3.5 kpc (figure [Tn|. Figure [THl reveals several dips and 
peaks in Vrot/o"z during the formation of the object indicat- 
ing several misaligned populations, which prevent the for- 
mation of a cold disk. The syste m bears some resem blance 
with polar bulge galaxies (see e.g. ICorsini et al.ll2012l ). 

Models starting in major orientations and developing 
strong misalignments are represented by C-Maj, the sur- 
face density projections of which are shown in the bottom 
row of figure [5] Two stellar components at an angle ~ 60 
degrees are clearly visible. Consequently, the circularity dis- 
tribution in figure [TT] shows two distinct peaks at e ~ 0.85 
and ~ 0.25. The latter represents the ring which forms at 
late times (see figure [TS]). In the gas surface density projec- 
tions, fiows are visible, which channel the material to the 
centre. The vertical profile at _R = 5 kpc (figure [TO} is flat 
for 1 2; I < 3.7 kpc before dropping exponentially, whereas the 
radial profile (figure [9]) resembles a double-exponential with 
a dominant central component at _R < 5 kpc and rings at 
-R > 10 kpc. These proper ties of the model revea l a similar- 
ity to polar ring galaxies (|Whitmore et al.lll990r ). 

5.5 How reorientation affects the shape of the 
potential 

In the context of the results already presented in this section, 
the question arises, how the shape of the final gravitational 
potential differs between models. Clearly, in central regions 
dominated by a disk, the potential is found to be oblate. 
Due to the interaction with the forming disk, the dark mat- 
ter halo_becoinesJess triaxi al, in agreement with the results 
of Kazant zidis et al.l (|201G| ). Moreover, the initial triaxial- 
ity of the halo, which in the dark matter only simulations 
shows little reorientation, is still imprinted at outer, halo- 
dominated radii. Consequently, for models such as AExt or 
CExtCosmo, the gas angular momentum of which was ini- 
tially aligned with the minor potential axis, the final disk 
axis is well aligned with the final minor axis of the poten- 
tial at all radii. For disk models, which end up in preferred 
disk orientation, but started with a different initial orien- 
tation, such as ARef-fMed or CExt-(-Arb, the situation is 
similar, yet the potential contours are less flattened, consis- 
tent with the fact that these disks are thicker. For models, 
where the flnal stellar angular momentum does not align 
with the preferred disk orientation, such as ARef-Med or 
CExt-Arb, there is also no alignment between the flnal po- 
tential ellipsoids at inner and outer radii. This is also true 
for compact disk models, that formed from gas which was 
initially aligned with the major axis. The disks, which are 
finally closest to an alignment with the major axis of the 
outer potential are AMaj and AMaj2C with an offset angle 
of ~ 30 degrees. 
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Figure 19. Orientation angle as a function of time for models ACosmo (top left)and CExtCosmo (top right) The initial gas sphere has 
been divided into 12 radial shells of i?gas/12 width each. For each shell the angle of the angular momentum vector ji(t) of all particles 
initially in shell i to the initial orientation is calculated at each output time t. The panels in the bottom row show the spin orientation 
in the same coordinate system for the gas content of 12 spatially fixed concentric shells of 10 kpc width for the A and C runs of CS09. 



5.6 Summary 

In conclusion, we have shown, that each of the two haloes 
we studied shows a preferred orientation for disks, which 
roughly agrees with the minor axis of the inner halo po - 
tential, in agreement with the results of lBailin et alj (|2005l ) 
for fully cosmological simulations. However a flip of the ini- 
tial angular momentum vector by 180 degrees produces a 
slightly worse (halo C) or unsuitable (halo A) orientation 
for a forming disk. This shows that it is not the shape of 
the potential alone that determines the preferred orienta- 
tion. Models with initial non-preferred angular momentum 
orientations tend to reorient to a preferred axis, producing 
thickened and flaring disks. In some cases they fail to settle 
to a stable orientation and do not end up with significant 
disks. There is no obvious simple criterion to explain this di- 
chotomy. Models starting with angular momentum parallel 
to the major halo axis also do not show a stable orientation. 
They form stable disks only if their angular momentum con- 
tent is low, so that they rapidly form compact objects, which 
cannot be destroyed by reorientation, which misaligns inner 
and outer components. Realistic thin disks with continu- 
ously forming populations are thus only possible for initial 
orientations close to the halo minor axis. Models with ini- 
tial orientations in between these cases are not capable of 



producing cold disks. Taking results of the previous sections 
into account, we have also shown that the orientations both 
of the central disk and of the outer gas can change by more 
than 90 degrees, despite the simple and coherent initial dis- 
tributions of our gas components. This is also true for the 
angle between the components. 



6 COMPARISON TO COSMOLOGICAL 
SIMULATIONS 

Until now, we have shown that, with our method, we are 
capable of introducing disks into Aquarius haloes A and C. 
Using the cosmological orientations of the disks a,t z — 1.3 in 
the simulations of CS09, we showed that for halo C the disk 
is in the 'right place', meaning that our model CExtCosmo, 
which starts with this orientation, produces a nice disk at 
z — 0. In halo A, the cosmological disk is in the 'wrong 
place', i.e. our model ACosmo, starting in the corresponding 
orientation, fails to produce a nice disk. We have identified 
reorientation of the angular momentum of infalling material 
and the angular momentum losses and misalignments con- 
nected to this phenomenon as the key factors in the failure 
of many of our models to build substantial disks. CS09 also 
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Figure 20. The angle between the angular momentum vector of all stars within i? < 10 kpc and the angular momentum vector of all 
dark matter within R as a function of R at z = 1.3 and z = 0.0 for models ARef (top left), CExtCosmo (top right) and the results of 
CS09 for haloes A (bottom left) and C (bottom right). The effect of varying the radius that defines the galaxy for halo A is indicated in 
the bottom left panel, where the blue line is for a 30 kpc galaxy radius. For all other models this radius definition plays a minor role. 



identified misaligned infall of gas as the problem destroying 
the existing disk structure in their resimulation of halo A. 

In this section we would therefore like to analyze the 
reorientation of gas in the runs of CS09 and compare them 
to our own. In these cosmological simulations the inflow of 
gas is continuous and more complex than in ours, it is not 
desirable to define fixed gas populations by their radial po- 
sition at a certain time. It is more instructive to follow how 
the angular momentum of gas evolves in fixed radial shells 
of 10 kpc width as is shown in figure [T91 where we compare 
the reorientation in models ACosmo and CExtCosmo to the 
reorientation in models A and C from CS09. 

ACosmo, as depicted in the top left panel of figure [T91 
shows strong and continuous reorientation for all baryons, 
gaseous and stellar. Reorientation is slightly delayed for the 
outer shells compared to the inner, but the resulting mis- 
alignments are small compared to the models discussed in 
the previous section. After t ~ 9 Gyr there is no overall 
detectable misalignment among the initially defined shells, 
although there is a small counter-rotating component as 
shown in figure [16] originating from temporary misalign- 
ments and continuous reorientation. As shown in figure [TH 
ACosmo is ~ 45 degrees offset from the orientation of most 
of the disk models in halo A at 2 = 0. During its evolution 



it evolves towards this orientation until f ~ 11 Gyr, when 
it is only ^ 20 degrees offset and then starts to turn away 
again. 

For the fully cosmological model we start the analysis 
at z = 2 (lower left panel). The reference orientation is the 
one of the stars within i? < 10 kpc at z — 1.3, which is the 
initial orientation of ACosmo. At 2 = 2 the orientation of 
the stellar angular momentum differs by less than 10 degrees. 
Even at z = 2 the gas out to _R = 120 kpc is not well aligned 
with the galaxy. Until 2; ~ 1 there is a clear misaligning 
process. From then on, gas in all shells shows an orientation 
~ 150 degrees offset from that of the galaxy at z = 1.3. The 
new orientation is also depicted in figure [HI It only differs 
by ~ 20 degrees from the initial orientation of model ARef, 
which explains why this orientation yielded the best results. 
The offset is caused by a stronger misalignment with the 
plane perpendicular to the major axis. 

The strong reorientation of the gas is also depicted in 
figure 10 of CS09, where the evolution of the angle between 
the stellar and the gas components within 27 comoving kpc 
is plotted. The fact that the components become aligned 
again at 2 = indicates that the inner stellar component 
reorients to the same direction, but with a delay. Moreover, 
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the central gas components (black lines in figure [T9)l are also 
offset fi-om the outer gas by up to ~ 50 degrees. 

The new orientation of the (outer) gas is also within 
~ 20 degrees of the orientation that model ACosmo shows 
at t '^ 10.5 Gyr, but ai z — the orientations differ by 
~ 70 degrees. This is due to the continual reorientation in 
model ACosmo. While the reorientation process is qualita- 
tively similar in these models, the final orientation of the 
galaxy is also dependent on the details of its assembly, as 
seen above for non-disk models with identical initial orien- 
tation but otherwise different parameters such as AReflSO 
and AExtlSO. 

For halo C the situation is less complicated. The ori- 
entation of the gas in the CS09 model, which was used as 
an initial condition for our model CExtCosmo at z = 1.3, 
behaves very similarly to CExtCosmo. It shows an overall 
reorientation by the same amount (~ 20 degrees) and com- 
paratively small misalignments, < 45 degrees. 

In figure [201 we compare the orientation of the galaxy 
(all stars within 10 kpc) in a given model to the orientation 
of the angular momentum of the dark matter (including sub- 
haloes) within i? as a function of radius R. Model ARef is 
initially set up to be aligned with the dark matter within 300 
kpc. As can be seen, this orientation is not, however, aligned 
with the dark matter within 200 kpc, but the stars and dark 
matter become aligned by z = where the alignment is, in 
fact, better than 30 degrees for all radii R < 500 kpc. 

In the fully cosmological halo A model the galaxy is 
aligned with the dark matter within R ~ 100 kpc ai z — 
1.3, but is strongly misaligned otherwise, especially at i? > 
150 kpc, where the angle reaches ~ 150 degrees. At z = 
0, there is no disk and thus the orientation of the galaxy 
is not well defined. Considering stars within 10 or 30 kpc 
yields different results as indicated by the black and blue 
lines. The stars at i? < 10 kpc are well-aligned with the 
central dark halo, whereas the stars at 10 < _R/kpc < 30 are 
not. Both populations are misaligned with the dark matter 
outside _R ~ 50 kpc, but the misalignment is significantly 
smaller for the stars at 10 < R/kpc < 30. 

For halo C the picture is different. Model CExtCosmo 
is initially misaligned with the dark matter angular momen- 
tum by an angle of ~ 50 degrees. This hardly changes by 
z = 0, when the inner 100 kpc have become more aligned, 
but the outer parts remain unchanged. In the fully cosmo- 
logical run the situation is strikingly different. In the inner 
~ 50 kpc, the dark matter both at z = 1.3 and at z — 0.0 
is strongly misaligned with the galaxy, almost 180 degrees 
near the centre. This explains why there is no clear correla- 
tion between dark matter and galaxy angular momenta in 
our models based on halo C. 

In summary, there are no large discrepancies between 
our models and fully cosmological galaxy formation simu- 
lations within the same haloes. In the absence of any de- 
structive events such as major mergers, a disk will survive 
over cosmological timescales, if it has formed in the preferred 
orientation for its halo. As we showed, this orientation does 
not depend on the potential alone (section 5), nor is there 
a clear correlation with the angular momentum of the halo 
(this section). The orientations of the halo p rincipal axes and 
the halo spin axis can change with time fsee lVera-Ciro et al.l 



I2OIII and lBett &: Frenkl2012l respectively) and so apparently, 
can the preferred disk orientation, as is clear from the halo 
A run of CS09, where a disk is destroyed by infiowing mat- 
ter oriented in a new preferred direction. Despite this we 
were able to find initial gas setups which lead to thin disks 
in this halo. For halo C there is no such orientation change 
and thus both our own and the fully cosmological models 
produce surviving disks. 



7 CONCLUSIONS 

We have presented a series of simulations of idealized SPH 
models for the formation and evolution of galactic disks 
within fuUy-cosmological ACDM haloes. At z — 1.3 we 
add rotating spheres of hot gas in approximate hydrostatic 
equilibri um to two dark mat ter haloes from the Aquarius 
Project (jSpringel et al.l l2008V A parametrized cooling law 
and standard prescriptions for star formation allow us to 
study the evolution of the combined baryon+CDM systems 
until z = 0. We study models with different orientations 
and amounts of baryonic initial angular momentum, differ- 
ent cooling timescales and different density profiles. This al- 
lows us to determine favorable and non-favorable conditions 
for the formation and survival of disks. 

Clearly, such simulations are not full models for the for- 
mation of disk galaxies. However, they allow us to study pro- 
cesses relevant to the formation of disks within triaxial and 
realistically evolving dark matter haloes, and hence explore 
the conditions under which massive disks can exist in ACDM 
haloes. This should help us understand the failure or suc- 
cess of various fully cosmological simulations of disk galaxy 
formatio n, which still suffer from a multitude of unc ertain- 
ties (e.g. ICovernato et al.ll2007l : ICuedes et al.l l 201ll . CS09, 



IScannapieco et al.ll2011h . A particular problem may be that 
these simulations tend to produce overly old stella r popula- 
tions (see the discussion ini Aumer fc Binnevll2009l ') and thus 
to under- predict the stellar mass formed since z '^ 1, the 
epoch generally thought to be best suited for disk formation 
due to the absence of major mergers at such late times. Our 
models are unable to reproduce full star formation histo- 
ries but yield interesting insights into galaxy formation at 
relatively recent times. 

Our setup has an inherent initial discrepancy between 
the triaxial, substructured and dynamically growing halo 
and the rotating quasi-equilibrium gas sphere. This leads to 
an initial phase, in which the gas adjusts to the potential. 
Thereafter it loses angular momentum to the dark matter 
as it cools to the centre and transforms the total central 
potential, and thus also the halo, into a more axisymmetric 
configuration capable of hosting a stable disk. We show that 
the shorter the timescale of this transformation, the more de- 
structive the effect on the forming disk, and thus the more 
dominant the resulting bulge and thick disk components. 
Haloes are expected to be n ear prolate after major merg- 
ers (e.g. iRomano-Diaz et al.ll20 09l before most of the stellar 
mass assembles into disks, so this process of bulge formation 
during halo transformation is not unrealistic. Romano-Diaz 
et al. also report the initial formation of an asymmetric, bar- 
unstable disk in their fuUy-cosmological simulations. Only if 
we align the initial angular momentum of the gas sphere 
with the major axis of the potential, can we suppress this 
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destructive process. However, the transformation from a pro- 
late potential to an oblate disk potential triggers misalign- 
ments between inner and outer baryonic components, which 
hinders disk survival. 

We have shown that our models are capable of pro- 
ducing, and thus that the haloes in consideration arc ca- 
pable of hosting, disks with realistic structural and kine- 
matic properties, which unlike in most cosmological simula- 
tions show realistic disk-to-bulgc ma ss ratios D/B > 4. As 
has been shown previously (Scannapi eco et al.ll2010l ). struc- 
tural and kinematic decompositions yield strongly differing 
results, with structural values exceeding kinematic ones by 
factors of a few. We show that our disks have rotation-to- 
dispersion ratios Kot/f^ monotonically decreasing with in- 
creasing age, which arises from continuous disk heating due 
to substructures and a monotonic decrease in birth veloc- 
ity dispersion connected to the transformation of the cen- 
tral potential. Due to limitations in resolution, our simula- 
tions are not capable of producing reali stic young thin ste l- 
lar disks with a^ < lOkms"^ (see also I House et al.ll201l[ ). 
but they agree qualitatively with age - velocity dispersion re- 
lations for the solar n eighborhood (|Holmberg et al.l l2009l : 
lAumer fc Binnevll2009l ). Our disks show truncations, which 
move outward with time and the radius of which increases 
with decreasing formation timescale and increasing angular 
momentum content. 

We find that the slower the formation of a model galaxy, 
the colder and more dominant is the final disk. This is the 
result of the combined effects of the suppression of bulge 
formation during potential sphericalization and of continu- 
ous disk heating being less efficient for younger populations. 
Bars develop in our model disks if their surface density is en- 
hanced by a higher total mass or lower total angular momen- 
tum. They strongly increase disk heating (see ISaha et al.l 
l20ld ). enhance the bulge fraction and prod uce prominent 
ring structures outside the bar regions (see lAthanassoulal 
1 19921 ) thus lowering the resulting Vrot/o"z by factors of a 
few. 

We show that the most stable, most extended and cold- 
est disks form from a gas angular momentum orientation 
that is aligned with the minor axis of the inner halo, as 
had been pr eviously found in cosmological simulations by 
ISailin et al.l (|2005V Out of the two such orientations given 
by the halo (180 degrees apart) one is clearly preferred. 
For one of our two haloes, the preferred direction is de- 
termined by the angular momentum of the outer halo in 
the initial conditions; by the end of the simulation all re- 
gions are (marginally) aligned with the central disk. The 
counter-rotating orientation is strongly disfavored and does 
not produce stable disks. For the second halo, there is no 
alignment between halo angular momentum and shape. The 
non-preferred minor axis orientation of the gas angular mo- 
mentum is the one, which makes the larger angle with the 
halo spin axis. Models with this initial orientation are able to 
host stable, yet thicker disks. In each halo we find a certain 
range of initial gas angular momentum orientations perpen- 
dicular to the halo major axis, for which a reorientation of 
the baryonic spin vector to one of the stable orientations oc- 
curs with the disk thickening but surviving. Models, where 
the gas angular momentum is initially aligned with the halo 
major axis can produce disks, but their orientation is not 
stable. Intermediate models cannot produce stable disks. Re- 



orientation especially heats the outer, lower surface density 
parts of disks, where self-gravity is unable to keep the stel- 
lar orbits aligned. Compact, high surface density disks are 
significantly less heated by orientation changes. 

As reorient a tion is common both for halo shapes 
( Vera-Ciro et al.l l201lh and for halo angular momenta 
teett fc FrenkI 120121 ') ■ the preferred orientation will evolve 
with time in ACDM haloes. This is directly connected to the 
well-known phenomenon of misaligned m atter infall during 
the a ssembly of haloes and galaxies (e.g. lOuinn fc BinnevI 
Il992r ). However, in our models, the initial gas distribution is 
spherically symmetric and its angular momentum is aligned 
at all radii. Despite this, the orientations of outer and inner 
gas and of stellar disks and also the angle between these 
components change by more than 90 degrees in many cases. 
Misalignment is thus a consequence of strongly nonlinear 
dynamical interactions in the final stages of galaxy assem- 
bly, not just of evolving asymmetries in the cosmological 
context, which feeds halo growth. 

Misaligned components ar e comm o n in o bserved galax- 
ies in the form of warps (jSancisil 19761) . polar ring s 
(|Whitmore et al.lll99Gl '). polar bulges jcorsini et al.l 1201^ 
or counter-rotating c omponents in early-type galaxies 
(jKrainovic et al.l 120081 '). Misaligned cosmological infall has 
been identified as a cause for a variety of phenomena, such 
as the de struction of disks (CS09), the formation of po- 
lar d isks ([Brook et al.l 120081 ) and the excitation of warps 
(Ros kar et al.ll20ld ). Interestingly, all of these phenomena 
also occur in our models. We qualitatively reproduce the 
misaligned infall found in the cosmological resimulation of 
halo A by CS09 starting from a spherical and coherently 
rotating gas distribution. The torques acting in our simu- 
lation must originate in the nonlinear structure of the dark 
h alo, suggesting a mec hanism different from the conclusions 
of lRoskar et al.l l|2010r ) who found negligible infiuence of the 
dark halo torques in comparison to anisotropic cosmological 
infall. 

All these processes of misalignment are associated with 
the loss of angular momentum (see also lRoskar et al.ll2010r ) 
and produce more compact, hotter objects. Thin extended 
disks can thus only form in stable orientations and thus 
in haloes for which the preferred disk orientation under- 
goes lit tle tempora l evolu tion. This conclusion is similar to 
that of I Sales et al.l (|201ll '). who found that disks form from 
the continuous accretion of material with coherently aligned 
angular momentum orientation and that misaligned accre- 
tion of gas tends to produce spheroids. The question arises 
whether ACDM offers enough haloes which fulfill these con- 
ditions to explai n the observed ab undance of disk galaxies 
in the universe. I Sales et al.l (12011! ) find a significant frac- 
tion of such disk-dominated systems in their model sample 
of galaxies, but the limitations of their simulations do not 
allow a full comparison with the statistical properties of ob- 
serv ed samples . Sem i-analytic models of galaxy formation 
(e.g. IGuo et al.ll201lh are helpful for studying galaxy popu- 
lations, but they do not capture the role of misalignments 
and other dynamical details that, as we have shown, play a 
major role in galaxy evolution. It would be desirable to add 
proper treatment of these processes to the models, however, 
we are not able to offer straightforward prescriptions, as we 
are not capable of drawing general conclusions from the pre- 
sented study of two haloes. 
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Considering our work and recent improvements in 
cosmoloRical resim u lations of individ ual galaxies (e.g. 
iGuedes et all 1201 ll : iBrook et al.l I2012D . realistic ACDM 
haloes appear to be capable of hosting realistic disk galax- 
ies. However, more work is needed to understand whether 
the population of ACDM haloes agrees in detail with the 
observed population of galaxies. 
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